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ABSTRACT 


Hot-wire  anemometry  measurements  were  made  in  the  12-in. 
supersonic  tunnel  of  the  von  Karman  Gas  Dynamics  Facility  to  deter¬ 
mine  the  level  of  flow  fluctuations  in  the  free  stream.  Supplementary 
acoustic  measurements  were  made  using  a  microphone  mounted  flush 
with  the  top  surface  of  a  Hat  plate.  Data  were  obtained  at  Mach  num- 
er  .  eynolds  number,  based  on  the  wire  diameter  (d  =  0.  0001  in.  ) 
ranged  from  5  0  to  24.0.  The  hot-wire  anemometry  results  presented 
are  based  on  the  mode  concept  of  Kovasznay  and  the  techniques  o' 
Morkovin  for  hot-wire  data  reduction.  The  free-stream  results,  which 
are  similar  to  those  published  by  Dr.  John  Laufer.  can  be  explained  by 
the  assumption  of  a  sound  source  traveling  in  the  tunnel  wall  boundary 
layer  at  a  supersonic  velocity  relative  to  the  free  stream.  The  accustic 
measurements  on  the  surface  of  the  flat  plate  revealed  pressure  fluc¬ 
tuation  levels  larger  than  those  of  the  free  stream  by  a  factor  of  at 
least  20. 
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SECTION  I 
INTRODUCTION 


The  fluctuating  variables  in  the  free  stream  of  a  supersonic  wind 
tunnel  test  section  have  been  accorded  much  interest  by  investigators 
of  such  allied  phenomena  as  boundary-layer  transition,  flow  separation, 
boundary-layer  control,  heat  transfer,  etc.  However,  few  measure¬ 
ments  of  supersonic  flow  fluctuations  have  been  published.  As  con¬ 
sideration  is  given  to  the  importance  of  high-speed  turbulence  and 
related  factors  in  the  design  and  performance  of  new  vehicles,  the 
capability  of  measuring  flow  fluctuations  may  become  a  requirement 
of  wind  tunnel  test  facilities.  The  capability  of  determining  flow  fluc¬ 
tuations  in  a  wind  tunnel  is  desirable  from  two  points  of  view:  the 
intrinsic  relation  of  the  fluctuating  quantities  to  phenomena  investi¬ 
gated  in  wind  tunnel  testing  and  the  inherent  role  of  the  fluctuations  in 
the  correlation  of  phenomena  among  various  wind  tunnels. 

Kovasznay  (Ref.  1)  has  shown  that  the  random  fluctuations  in  a 
supersonic  flow  are  composed  of  three  types  or  modes  of  disturbance: 
the  vorticity,  entropy,  and  sound-wave  modes.  Morkovin  ( Ref.  2)  has 
discussed  the  possible  sources  of  the  different  modes.  Briefly, 
Morkovin  has  reasoned  that  the  vorticity  and  entropy  fluctuations  in 
the  free  stream  have  their  origins  in  the  tunnel  settling  chamber  or 
upstream,  and  that  the  sound  fluctuations  are  generated  within  the 
boundary  layers  on  the  tunnel  walls.  Morkovin  points  out  that  vorticity 
(turbulent  velocity)  fluctuations  are  effectively  suppressed  by  a  large 
speed  ratio  across  the  nozzle.  He  also  comments  that  the  decay  of 
entropy  fluctuations  (temperature  spottiness)  is  accelerated  by  mixing 
upstream  of  the  settling  chamber  and  by  screens  in  the  settling  cham¬ 
ber.  It  is  concluded  that  the  magnitude  of  free-stream  vorticity  and 
entropy  fluctuations  in  the  supersonic  test  section  are  to  a  large  extent 
subject  to  control.  On  the  other  hand,  it  does  not  appear  that  design 
considerations  can  greatly  alter  the  level  of  the  sound  mode  of  fluctua¬ 
tions  in  the  free  stream. 

Some  years  ago  (in  1954)  the  importance  of  pressure  fluctuations 
(sound)  emanating  from  turbulent  boundary  layers  on  the  walls  of 
supersonic  wind  tunnels  was  suggested  by  Laufer  (Ref.  3).  Subsequent 
work  by  Laufer  (Refs.  4  and  5),  using  hot-wire  anemometry  techniques, 
established  the  level  of  sound  fluctuations  in  the  free  stream  of  one 
wind  tunnel.  More  recently  the  influence  of  aerodynamic  noise  radiated 
from  tunnel  wall  boundary  layers  upon  transition  of  the  bcjundary  layer 
of  test  models  has  been  investigated  by  Pate  and  Schueler  (Ref.  6). 
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tlla  pre!?nt  rep°rt  presents  results  of  measurements  of  the  fluc¬ 
tuating  variables  in  the  free -stream  flow  field  of  one  of  the  supersonic 

used  innpatS  ThT6!!  °i  °f  ^  V°n  Karman  Gas  Dynamics  Facility  (VKF) 

ments  it  has  he  1Sc!,“eleIf.  *  StUdieS’  In  conjutiction  with  these  measure- 
ments,  it  has  been  the  objective  of  this  investigation  to  develop  a  caoa- 

bihty  within  the  VKF  to  use  the  hot-wire  anemometer  as  a  measurement 
tool  in  high-speed  flows.  Although  the  heated  fine  wire  is  not  a  con- 

roS  tfiTs6  t0nUSe  andfthK  aSSOciated  data  analysis  cannot  be  called 
’  .  he  8“f1.1  Slze  of  the  measuring  element  and  the  ability  of  the 
wire  with  its  associated  electrical  equipment  to  respond  to  high 
frequency  fluctuations  make  the  hot-wire  anemometer  the  principal 
.nsttM  now  in  use  studying  high-speed  no.  ttuotuatioi  ' 

Hot-wire  measurements  were  made  at  Mach  number  4  at  three  unit 
S  numbers-'  °-  °-  «•  »"d  «•  *  million  per  inch.  Supp!2 

na'  pufe  W"e  mad'  “slng  *  microphone  mounted  on  a 


SECTION  II 
APPARATUS 


2.1  WIND  TUNNEL 

Tunnel  D  is  an  intermittent,  variable  density  wind  tunnel  with  a 
manually  adjusted,  flexible  plate-type  nozzle  and  a  12-  by  12-in  test 
section.  The  tunnel  can  be  operated  at  Mach  numbers  from  1.  5 ’to  5  at 
stagnation  pressures  from  about  3  to  60  psia  and  at  an  average  s  agna 
Uon  temperature  of  about  70»F.  A  description  of  the  tunnel  and  airflow 

tCh?T?nnerDni0rmat/0n  may  ^  fOUnd  in  Ref>  7‘  An  “Oration  showing 
he  Tunnel  D  geometry  is  presented  in  Fig.  1  (Appendix  I)  It  is  of  g 

interest  to  note  that  Tunnel  D  is  a  modified  copy  of  the  12-in  super¬ 
sonic  wind  tunnel  which  was  previously  in  operation  at  the  Jet  Propul¬ 
sion  Laboratory,  California  Institute  of  Technology.  P 

nn  h.0t',wire  measurements  of  the  present  investigation  were  made 

on  the  centerline  of  the  tunnel  at  the  centerline  of  the  viewing  windows 
or  approximately  55.  1  in.  downstream  of  the  throat.  At  thisf  location 

1^  in  Wc^UdnedoarydlayerSt0/  *?*  tUnnel  WallS  is  aPP™*imately  1.  0  to 
(Ref  7)  dePendl"g  ‘he  free-stream  unit  Reynolds  number 


nozzfeThfoaf  thS^'th  ’  ““  COntraction  ratio  upstream  of  the  Tunnel  D 
,.  .  u  that  1S  the  cr°ss -sectional  area  of  the  stilling  chamber 
lvided  by  the  area  of  the  throat,  is  approximately  84  to  1.  S 


2 


AE  DC-TR-71  -143 


2.2  HOT-WIRE  PROBES 

The  hot-wire  probes  used  in  this  investigation  (Fig,  2)  were  fabri¬ 
cated  by  the  VKF  and  were  similar  in  design  to  probes  used  successfully 
in  supersonic  flow  by  Laufer  and  Vrebalovich  (Kefs.  4,  8,  and  9). 
Needles  were  glued  to  a  razor  blade  section  attached  to  a  mounting  strut. 
The  needles  extended  0.  04  in.  upstream  of  the  leading  edge  of  the  blade, 
with  a  nominal  spacing  of  0,  015  in.  at  the  tips.  Probes  of  a  different 
design,  with  needles  extending  0.  7  in.  ,  which  were  used  in  earlier  VKF 
hot-wire  measurements  in  low  density  separated  flow,  were  tried 
initially  in  the  present  case,  but  the  wire  supports  were  not  suffici¬ 
ently  stable  relative  to  one  another  so  that  the  wires  on  these  probes 
failed  upon  exposure  to  the  free  stream. 

Platinum- 10  percent  rhodium  wire  drawn  by  the  Wollaston  process 
was  used  for  the  sensing  element  of  the  hot-wire  probes.  The  silver 
sleeve  of  the  Wollaston  process  which  surrounded  the  fine  wire  was 
removed  by  etching  in  nitric  acid  using  a  small  electric  de-plating  cur¬ 
rent.  After*  removing  the  sleeve,  the  wire  was  soft-soldered  to  the 
probe  needles,  leaving  slack  to  reduce  tension  in  the  wire  under  load¬ 
ing  (see  Fig.  2c).  Each  wire  was  annealed  by  heating  it  in  still  air  for 
approximately  five  minutes. 

The  manufacturer's  nominal  value  for  the  wire  diameter 
(d  =  1,  0  x  10  4  in.  )  was  accepted  for  the  data  reduction.  According  to 
the  manufacturer,  the  wire  was  drawn  to  a  resistance  tolerance  of 
±5  percent,  not  to  a  diameter  measurement  tolerance.  A  value  for  the 
resistivity  of  the  annealed  platinum- 10  percent  rhodium  wire  (111  ohms 
pei*  circular  mil  foot)  was  obtained  from  the  manufacturer  and  was  used 
to  calculate  the  length  of  each  wire  from  the  measured  resistance.  The 
two  probes  used  for  the  measurements  had  wire  aspect  ratios  {£/ d)  of 
approximately  140  and  160. 

The  relationship  between  the  resistance  and  temperature  of  a  heated 
fine  wire  may  be  stated  by  the  equation: 

«w  =  «f[>  -  ♦  0,<VTf)*] 

where  the  subscript  w  refers  to  the  heated  conditions  and  the  subscript  f 
to  the  chosen  reference  conditions.  A  sample  of  the  wire  used  in  the 
present  measurements  was  check-calibrated  against  the  nominal  value 
for  obtained  from  Morkovin  (Ref.  10)  (<*f  =  8.  9  x  10‘4  per  °R  for 
Tf  =  49  2°R ),  The  nominal  value  was  acceptable  within  the  estimated 
uncertainties  of  the  calibration.  The  value  for  j3f  obtained  from  Morkovin 
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(nominal  P{  =  -4. 8  x  10  8  per  °R2)  was  of  negligible  significance  for  the 
range  of  applicable  temperatures  encountered  in  the  present  study. 


2.3  MICROPHONE 

Supplementary  pressure  fluctuation  measurements  were  made  using 
Bruel  and  Kjaer.  Model  No.  4136,  0.  25 -in.  -diam  condenser  microphone 
cartridge  mounted  in  a  flat  plate,  as  shown  in  Fig.  3.  The  microphone 
system  consisted  of  the  following  Bruel  and  Kjaer  components  connected 
to  the  cartridge:  a  flexible  adapter,  type  UA0122;  a  cathode  follower 
type  2615;  and  the  power  supply,  type  2801.  The  system  had  a  flat 
,±2  d' eclbel  (db))  frequency  response  at  atmospheric  pressure  from 
30  Hz  to  70  kHz  and  a  dynamic  pressure  range  of  from  70  to  180  db 
referenced  to  a  sound  pressure  of  0.  0002  microbar.  The  sensitivity  of 
the  microphone  system  in  terms  of  output  voltage  to  applied  pressures 
was  determined  by  VKF  calibration  to  be  0.  070  millivolt /microbar. 

The  flat -plate  model  was  designed  to  accommodate  a  second  micro¬ 
phone  mounted  internally  (Fig.  3)  to  monitor  the  effects  of  model  vibra¬ 
tion  on  microphone  response.  Inasmuch  as  these  effects  had  previously 
been  found  negligible  (Ref.  6)  under  similar  test  conditions,  the  second 
microphone  was  omitted  from  the  present  study.  The  air  mass  within 
the  condenser  microphone  provides  mechanical  damping  for  the  micro¬ 
phone  diaphragm.  Design  considerations  make  the  damping  critical  at 
an  ambient  pressure  of  one  atmosphere,  in  order  to  give  the  widest 
frequency  range  of  flat  response  to  sound  pressure  fluctuation  at  that 
condition.  When  the  microphone  is  operated  in  an  environment  of  sub- 
atmosphenc  pressure,  the  frequency  response  is  altered,  especially  at 
the  higher  frequencies,  by  a  resonance  peak,  and  the  output  signals 
which  are  centered  about  the  resonance  frequency  (45  kHz  in  the  present 
case)  introduce  errors  into  the  sound  pressure  measurements.  For 
present  purposes,  however,  the  resonance  contribution  was  determined 
by  playback  of  the  tape-recorded  microphone  output  through  a  Spencer- 
Kennedy  Model  302  variable  electronic  bandpass  filter  set  to  pass 
signals  from  45  to  55  kHz.  (The  attenuation  at  the  filter  cutoff  fre¬ 
quencies  was  1.  8  db  with  a  rejection  rate  of  18  db  per  octave,  accord¬ 
ing  to  the  manufacturer's  specifications. )  The  mean -square  of  the 
resonance  contribution  was  read,  using  the  true  rms  meter,  for  each 
measurement,  and  the  result  was  then  subtracted  from  the  mean-square 
ot  the  total  microphone  output. 
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2.4  HOT  WIRE  AND  MICROPHONE  INSTRUMENTATION 


The  hot-wire  anemometer  measurements  were  made  using  com¬ 
mercially  available  constant-current  equipment  (Fig.  4)  of  Shapiro  and 
Edwards  design  (Model  50). 

The  amplifier  (Model  A-50C)  has  a  design  frequency  range  of  from 
1  Hz  to  320  kHz.  Coupling  circuits  are  provided  within  the  amplifier 
to  permit  control  of  the  low  frequency  and  high  frequency  cut-off  char¬ 
acteristics.  The  amplifier  has  a  maximum  gain  of  50,  000  between  in¬ 
put  and  output,  with  the  compensation  network  off.  Maximum  gain  when 
fully  compensated  is  2.  5  x  10°.  The  gain  of  the  amplifier  is  variable 

a  range  of  1028:1  in  21  steps  of  0.  707:1  ±  1  percent  (approximately 
3  db)  each.  The  maximum  uncompensated  gain  was  checked  at  least 
once  each  day  of  operation. 

A  compen  mating  network  is  incorporated  in  the  amplifier  to  offset 
the  effect  of  the  thermal  inertia  of  the  hot-wire  at  higher  frequencies. 
The  network  is  capable  of  compensating  wires  with  time  constants 
between  0.  05  and  22.  0  msec.  The  ratio  between  basic  amplification 
and  the  maximum  amplification  of  the  compensating  circuit  is  500  for 
any  setting  of  the  time  constant. 

W.re-heating  currents  can  be  adjusted  over  the  range  of  from  1 
to  300  ma.  The  current  was  measured  using  a  precision  potentiometer 
with  ranges  of  1.  0,  0.  1,  and  0.  01  amp  full  scale  with  a  precision  of 
±0.  1  perce:  t  of  full  scale. 

The  resistance  of  the  wire  was  measured  using  a  precision  bridge 
with  a  precision  of  ±0.  01  ohms  or  ±0.  05  percent  of  the  reading,  which¬ 
ever  is  greater. 

The  ratio  of  signal-to-noise  was  improved  by  using  an  input  trans¬ 
former  which  increased  the  input  signal  voltage  to  the  amplifier  by  a 
factor  of  five.  The  transformer  has  two  frequency  ranges:  25  Hz  to 
40  kHz  and  1  kHz  to  320  kHz.  In  the  present  investigation,  data  were 
obtained  for  both  ranges. 

The  rms  voltage  fluctuations  of  the  wire  were  fed  from  the  ampli¬ 
fier  to  a  Briiel  and  Kjaer  true-rms  meter,  type  2409.  Each  reading 
(which  had  an  uncertainty  of  ±2.0  percent)  was  corrected  for  amplifier 
electronic  noise  determined  for  the  individual  control  settings  of  the 
amplifier.  Microphone  output  was  measured  using  the  same  meter. 
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Magnetic  tape  recordings  of  the  hot-wire  anemometer  amplifier 
output  were  made  using  an  Ampex  Model  FR  1300  analog  tape  recorder. 
The  hot-wire  signals  were  recorded  simultaneously  on  a  direct  channel 
and  an  FM  channel  having  frequency  responses  from  300  Hz  to  300  kHz 
and  dc  to  20  kHz,  respectively.  Signals  from  the  microphone  on  the 
flat  plate  were  recorded  for  tunnel  flow  conditions  which  duplicated  (on 
separate  runs)  those  for  the  hot-wire  data.  Microphone  signals  were 
recorded  only  on  a  direct  channel. 

Power  spectra  of  the  hot-wire  response  signals  and  of  the  micro¬ 
phone  signals  were  determined  using  a  Honeywell  9300  automatic  wave 
analyzer  system.  This  system  is  designed  to  analyze  signals  in  the 
frequency  range  from  2  Hz  to  40  kHz,  Available  analysis  bandwidths 
are  1,  5,  10,  and  50  Hz.  The  analysis  of  signals  of  higher -frequency 
content  may  be  obtained  through  tape  playback  speed  reduction.  In  the 
present  case,  selected  hot-wire  signals  were  analyzed  to  give  power 
spectra  over  the  frequency  range  from  20  Hz  to  120  kHz. 


SECTION  III 

DISCUSSION  OF  HOT  WIRE  MEASUREMENTS 


3.1  DATA  ACQUISITION  AND  REDUCTION  CONSIDERATIONS 

The  analysis  of  the  structure  of  turbulence  in  a  compressible, 
viscous,  heat -conductive  gas  by  perturbation  techniques  has  been  re¬ 
ported  by  Kovasznay  (Ref.  1)  and  Chu  and  Kovasznay  (Ref,  11), 
Kovasznay 's  linearized  analysis  (Ref.  1)  identified  three  types  or  modes 
of  fluctuations:  vorticity,  entropy,  and  sound.  In  that  analysis  the 
three  modes  are  treated  as  completely  decoupled;  that  is,  energy  trans¬ 
fer  from  one  mode  to  another  is  not  allowed.  Chu  and  Kovasznay 
(Ref.  11)  generalized  the  analysis  to  include  interactions  between  modes 
in  the  form  of  bilateral  (second-order)  terms.  Fortunately,  in  a  small 
space-time  domain  without  solid  boundaries,  such  as  the  close  neighbor¬ 
hood  of  a  fine  heated  wire,  the  interactions  are  of  a  weak  nature,  and 
therefore,  may  be  neglected.  The  concepts  of  Refs.  1  and  11  serve  as 
a  basis  for  the  interpretation  of  the  fluctuation  measurements  made 
using  the  hot-wire  anemometer. 

The  basic  relation  for  the  response  of  the  hot  wire,  immersed  in  a 
flow  field,  and  the  associated  electrical  system  to  fluctuations  of  the 
flow  is  stated  as  follows: 
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where^Ae  is  the  instantaneous  voltage  change  from  the  mean  voltage 
drop  e  across  the  wire;  ft  is  the  time  constant  of  the  wire;  t  is  time; 

Ap»  Au,  AT0  are  the  fluctuating  components  of  the  mean  flow  den¬ 
sity  p,  velocity  u,  and  total  (stagnation)  temperature  T0>  respectively; 
and  Aep,  Aeu,  and  AeT  are  the  sensitivity  coefficients  of  the  heated 
wire.  The  time  constant  and  the  sensitivity  coefficients  are  dependent 
on  the  properties  of  the  wire,  of  the  electrical  system,  and  of  the  mean 
flow.  The  second  term  on  the  left  side  of  the  equation  is  effectively 
made  to  vanish  by  proper  compensation,  within  the  electrical  system, 
for  the  thermal  inertia  of  the  wire. 

In  order  to  infer  flow  fluctuation  quantities,  the  sensitivity  coeffi¬ 
cients  must  be  evaluated  either  through  a  calibration  technique  where 
flow  conditions  are  known  or  by  relating  the  heat -transfer  character¬ 
istics  of  the  wire  to  the  electrical  circuit  parameters  through  empirical 
correlations.  The  latter  technique  is  most  often  used  because  of  the 
short  life-expectancy  of  a  given  wire.  Consideration  will  now  be  given 
to  the  derivation  of  expressions  for  the  sensitivity  coefficients. 


The  energy  balance  for  the  heated  wire  without  end  losses,  aver¬ 
aged  over  the  length  of  the  wire,  is  simply  stated  by  the  equation 

dK.'dt  =  W  -  II 

where  E  is  the  heat  stored  in  the  wire,  W  represents  the  thermal  energy 
generated  in  the  wire  per  unit  time,  and  H  represents  the  rate  of  heat 
loss  to  the  fluid  flowing  around  the  wire.  The  heat  stored  in  a  wire, 

E  =  CT  V  at  temPerature  Tw  is  viewed  as  constant  when  the  wire  is 

properly  compensated  for  thermal  lag  resulting  from  the  finite  heat 
capacity  of  the  wire.  The  power 

W  --  J 1 2  H w  =  JI2  Rfll  +  cqO’w-Tf)  +  Pf(Tw-Tf)2] 

is  generated  by  the  heating  current  I,  but  W  also  ’epends  on  the  tem¬ 
perature  Tw  of  the  wire  relative  to  a  reference  temperature  Tf  and  on 
the  calibration  of  the  wire  which  yields  the  constants  %  cvf,  and  0f 
(J  is  a  units  conversion  factor).  The  heat  loss  H  =  Nu0  irS.  k0(Tw  -  Tr) 
is  stated  in  terms  of  the  Nusselt  number  Nu0  =  h  d/k0  and  the  "equilib¬ 
rium"  or  recovery  temperature  Tr  =  r]  T0  of  the  unheated  wire. 

In  order  to  relate  the  heat  loss  of  a  hot  wire  to  the  conditions  of 
the  flow,  it  is  necessary  to  know  the  steady-state  heat  loss  relationship 
applicable  to  the  given  flow  field.  Such  relationships  are  usually  stated 
in  terms  of  the  Nusselt  number  (Nu0)  as  a  function  of  Mach  number  (Mm) 
and  Reynolds  number  (Re0  =  d/n0)  where  the  air  thermal 
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conductivity  kQ  and  viscosity  /i0  are  based  on  the  stagnation  tempera¬ 
ture.  *  In  supersonic  flow,  this  reference  temperature  has  the  advantages 
of  remaining  unchanged  across  shock  waves  and  of  being  a  good  approxi¬ 
mation  for  the  free-stream  temperature  (T2)  behind  a  normal  shock  wave 
(or  an  almost-normal  wave,  such  as  the  shock  wave  immediately  up¬ 
stream  of  the  blunt  wire).  This  fact  becomes  particularly  advantageous 
when  it  is  realized  (Ref.  12)  that  flow  conditions  behind  the  detached 
shock  wave  and  in  the  vicinity  of  the  wire  are  only  weakly  dependent  on 
the  free-stream  Mach  number.  Further,  when  calculating  the  energy 
loss  from  the  wire,  it  is  necessary  to  know  the  functional  dependence  of 
the  temperature  recovery  factor  (17)  on  the  relevant  flow  parameters. 

The  required  relations  can  be  obtained  by  calibrating  the  wire  in 
the  Hew  field  to  be  investigated.  It  should  be  noted  that  this  calibration 
is  concerned  with  mean  values  of  the  variables,  not  with  fluctuations. 
Alternatively,  quantitative  results  can  be  obtained  at  the  sacrifice  of 
some  accuracy  by  reliance  on  existing  correlation  functions  (Refs.  9 
and  10).  Relationships  have  been  established  for  infinite -length  wires 
as  a  result  of  the  work  of  many  investigators  and  have  been  formulated 
empirically  for  wide  ranges  of  Mach  number  and  Reynolds  number  by 
Dewey  (Ref.  13)  and  others.  In  the  present  investigations  empirical 
relations  were  used  in  lieu  of  direct  calibration  of  each  wire  because  of 
the  short  life-expectancy  of  a  fine  wire  in  supersonic  flow.  For  sim¬ 
plicity,  tiie  data  of  the  present  investigation  have  been  reduced  using 
information  presented  by  Vrebalovich  (Ref.  9),  who  has  expressed 
Nusselt  number  for  (infinite -length)  wires  in  supersonic  flow  (M  >  1.2) 
as  a  function  of  Reynolds  number  only  and  has  presented  a  table  of 
experimental  data  relating  the  recovery  temperature  factor  to  Reynolds 
number  only  (see  Fig.  5). 

Thermal  lag  due  to  the  finite  heat  capacity  of  a  wire  causes  the 
response  of  the  heated  wire  to  flow  fluctuations  to  decrease  at  higher 
frequencies.  The  compensation  network  within  the  anemometer  ampli¬ 
fier  circuit  permits  electrical  compensation  for  the  thermal  inertia  of 
the  heated  wire;  that  is,  the  response  of  the  amplifier  is  made  to  in¬ 
crease  at  higher  frequencies.  As  a  consequence,  however,  electronic 


In  the  general  case  of  heat  transfer  from  an  infinite -length  wire, 
the  Nusselt  number  is  a  function  of  Reynolds  number,  Mach  number, 
and  temperature  loading  (Tr)  for  a  given  gas,  with  constant  Prandtl  num¬ 
ber  and  constant  ratio  of  specific  heats  (7).  The  experiments  of  Laufer 
and  McClellan  (Ref.  12)  have  shown  that  the  temperature  loading  is  a 
weak  parameter  and  that  the  slope  of  the  Nusselt  number  versus  Reynolds 
number  correlation  varies  only  slightly  over  a  wide  range  of  temper  ature 
loadings . 
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noise  is  also  amplified,  and  the  signal -to-noise  ratio  decreases.  To 
minimize  the  thermal  lag  and  obtain  a  useful  signal-to-noise  ratio,  the 
diameter  and  length  of  the  wire  must  be  quite  small.  It  must  be  noted, 
however,  that  in  order  to  obtain  correct  high  frequency  response  it  is 
necessary  to  minimize  temperature  gradients  along  the  wire  by  re¬ 
quiring  a  large  length-to-diameter  ratio.  The  availability  of  an  input 
transformer  to  boost  the  strength  of  the  signal  from  the  wire  was  found 
to  be  valuable  in  obtaining  reasonable  signal-to-noise  ratios  for  the 
supersonic  flow  field  of  the  present  study. 


The  time  constant  for  thermal  lag  compensation  was  determined 
for  each  data  point  by  the  square-wave  technique  (Ref.  14).  Some  un¬ 
certainty  is  associated  with  this  method  of  experimental  time  constant 
determination  because  of  the  superposed  turbulence  signal  and  amplifier 
noise.  The  computational  schemes  of  Morkovin  and  Phinney  (Ref.  15) 
and  Vrebalovich  (Kef.  9)  were  used  to  minimize  the  uncertainties. 
Morkovin  and  Phinney  (Ref.  15)  show  that  the  ratio  of  the  time  constant 
of  the  correct  compensation  setting  to  that  of  an  incorrect  setting  is 
equal  to  the  ratio  of  the  correct  value  of  the  rms  voltage  output  to  the 
incorrect  value,  provided  that  the  energy  density  contribution  to  the 
overall  signal  of  the  uncompensated  wire  is  small.  The  difficulty  with 
this  correction  procedure  lies  in  determining  the  "correct"  time  con¬ 
stant  value.  To  overcome  this  difficulty,  Vrebalovich  suggests  (Kef.  9) 
calculating  the  heat  capacity  (Op)  of  the  wire  for  each  point  in  a  set  of 
data,  obtaining  the  average  heat  capacity  (C'p)  for  the  set,  and  assuming 
the  average  value  to  be  the  "correct”  value.  The  "correct"  time  con¬ 
stant  is  then  calculated  from  and  applied  in  the  correction  technique 
of  Morkovin  and  Phinney  (Ref.  15).  Discussion  of  the  effect  of  this  cor¬ 
rection  will  be  deferred  until  a  later  section  of  this  report. 

Inasmuch  as  the  electric  circi  :t  which  supplies  current  to  heat  the 
wire  of  the  ’constant-current"  anemometer  has  a  finite  impedance,  the 
current  is  not  truly  constant  when  the  wire  resistance  is  made  to  vary 
by  fluctuations  in  the  flow.  The  variations  in  the  voltage  (e)  and  thermal 
energy  (W)  are  small,  but  they  must  be  taken  into  account  using  a  finite - 
circuit  parameter  in  the  equation  for  the  conservation  of  heat. 

The  factors  discussed  above  which  are  related  to  the  heat  balance 
of  the  hot  wire  have  been  considered  by  Morkovin  (Ref.  10)  in  obtaining 
expressions  for  the  sensitivity  coefficients.  Assuming  proper  compen¬ 
sation  so  that  dE/dt  =  0  and  introducing  small  perturbations  into  the 
steady-state  energy  equation  (W  *  H)  in  logarithmic  form,  Morkovin 
derived  the  following  expressions: 
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where  E'  is  a  finite-circuit  parameter  and  where  A;  =  05OfnR  /3fnI) 
is  the  local  logarithmic  slope  of  the  resistance-current  variaUonand  Se*1 
’“*Z  *  h  f"f  ““  cr,r°“ed  h'*'1"8  »'  -hue  the  res,  of tte 

nfuch  «  Jhe  “1  g  he<  n°W  parameters-  are  maintained  constant.  Inas- 

field  a«  will  C»!  u  uCe  may  Vary  Wlth  Wire  location  m  the  flow- 
fieid  as  wen  as  with  the  heating  current,  the  slope  A'  must  be  deter- 

SS  S)  '“A-T.Thf  Cl,  lear‘  “  ,h=  description  by  Morkovip 
parameter  ,,  W  Uy  Hnearized'  locally  calibrated  overheat 

The  various  logarithmic  derivatives  which  appear  explicitly  in  the 
above  equations  for  the  sensitivity  coefficients  are  determined  from^he 

h  T  Nusselt  number  and  temperature  recovery  factor  obtained 
bv  calibration  or  through  use  of  the  correlations.  It  must  be  pointed  out 
here  that  a  significant  portion  of  the  heat  generated  in  the  hot  wire  is 

SSfaTrTl  ^TisVjoT^rre^ 

Z ‘^“T. P  ““  m',S-  9  “d  ■»>  «.  account  for  end  ZZ.il  to.ee, 
Tornm.rck  »•*  developed  by  Kov.een.y  and  ' 

Ref  T<  End  loj  again  presented  in  detail  by  Kovasznay  in 

been  ^  M  haVe 

inHpn*n^n  tht  wire  ft Maelt  number  and  temperature  recovery  factor  are 

e.Ztro,*rh  rf"Z  “  in  no.  field.,  .he  “re 

8itivity  coefficients  for  density  and  velocity  are  equal: 

Aep  .  Aeu  -  Aem 

EZfST  **  ‘he  8en8!tivity  t0  ma8S  now  fluctuations.  The  basic  equa- 
n  for  the  response  of  a  properly  compensated  hot  wire  , 


tp(¥)~  Ae“(?)+  Aer(x) 
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becomes,  then,  for  a  supersonic  flow  field, 
Ae  =  -Aein^~^  +  Act 
with  mass  flow  fluctuations. 


£  ■  (¥>  •  (t) 

and  total  temperature  fluctuations, 


The  mean-square  of  the  voltage  fluctuation  output  of  the  heated  wire 
is  given  by 


where  the  tilde  C)  indicates  an  rms  fluctuation  and 


(Am  HAT  0) 

R  I  -  - — 

in  I  f| 

is  the  correlation  coefficient.  Measurements  made  at  three  distinct 
wire  heating  currents  will  yield  three  equations  for  simultaneous  solu¬ 
tion  for  the  three  unknowns  (m,  T0,  and  RmT^  but  uncertainties  in  the 
measurements  are  likely  to  produce  scatter  in  the  results.  The  effects 
of  the  uncertainties  can  be  minimized  by  using  a  greater  number  of 
measurements  and  a  graphical  solution  for  the  unknowns.  This  manner 
of  presentation  of  hot-wire  data  was  introduced  by  Kovasznay  (Ref.  1) 
for  supersonic  flow-field  measurements  and  involves  only  the  rewriting 
of  the  basic  equation  in  the  form: 


and  the  determination  of  values  for  the  coefficients  which  give  the  best 
curve  among  the  measurements  plotted  in  the  form  of  (e/Aep)  versus 
(Aem/AeT).  Further  discussion  of  the  use  of  the  Kovasznay-type  of 
mode  diagram  will  be  found  in  Section  3.  2. 

While  the  wire  responds  to  fluctuations  of  the  density,  velocity, 
and  temperature,  the  independent  fluctuations  which  characterize  the 
turbulent  flow  field  are  the  vorticity  (turbulence),  entropy  (temperature 
spottiness),  and  sound  modes.  These  modes  are  related  to  the  density. 
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velocity,  and  total -temperature  sensitivity  coefficients  given  above  by 
the  following  expressions  (Ref.  11)  [the  expanded  equations  are  given 
in  Table  I  (Appendix  II)]: 


Aer  ~  -Aeu  +  6  Ae 
M  T 

Aeff  -  Aep  +  oM  Ae? 

°  X 

Ae^  -  -A t0 —  Ae„  +  a  (y-l)(l +nxM)Aej 

'  M  M 

where  AeT,  AeQj  and  Ae^  are  the  vorticity,  entropy,  and  sound  sensi¬ 
tivity  coefficients,  respectively,  and  nx  is  the  direction  cosine  of  the 
normal  to  the  sound  plane  wave  front  relative  to  the  flow  direction.  If 
the  flow  field  is  characterized  by  more  than  one  sound  direction,  that 
is.  if  two  or  more  significant  sound  sources  with  distinct  orientations 
exist  in  the  flow  field,  additional  sound  sensitivity  coefficients  v*ould  be 
required,  one  for  each  sound  wave  direction. 

The  hot-wire  voltage  fluctuation  equation  can  now  be  rewritten  as 
Ae  =  Ae^r)  +  A  ea(o)  +  Ae^) 


The  mean-square  voltage  fluctuation  output  of  the  hot  wire  is  given 

by 

(7)2  =  A^T)2  ♦  Ae,2*?)2  +  Ae„w  +  2AefAeffOO(<r)  Rro 
+  2Ae0rAen((7 ) (  fr )  +  2Ae^Ae^  ft )  (  r) 

where  (~)  again  denotes  the  rms  value  of  the  fluctuation  and  each  R  is 
a  correlation  coefficient  identified  by  its  subscripts.  In  principle,  six 
equations  in  six  unknowns  (the  three  fluctuations,  t,  a,  and  rf  and  the 
three  correlation  coefficients)  could  be  obtained  for  simultaneous  solu¬ 
tion  if  measurements  were  made  at  six  distinct  wire  heating  currents. 
(This  would  not  be  true  if  the  orientation  nx  of  the  sound  wave  were 
multivalued,  since  additional  sensitivity  coefficients  and  correlation 
coefficients  would  be  introduced.)  In  reality,  however,  the  determina¬ 
tion  of  the  sensitivity  coefficients  and  the  measurement  of  the  rms 
voltage  e  have  uncertainties  which  rule  out  this  method  of  solution 
(Ref.  10). 

Depending  on  the  nature  of  the  flow  in  which  the  measurements  are 
made,  usually  some  of  the  terms  in  the  equation  for  (  e  )2  can  be  ne¬ 
glected,  being  of  the  order  of  the  uncertainties  or  small  relative  to  the 
remaining  terms.  Specifically,  the  correlation  coefficients  (R^t  and 
RCT7r)  are  almost  always  negligible  because  the  nature  of  sound 
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propagation  differs  basically  from  the  manner  of  vorticity  and  entropy 
propagation  (Refs.  1  and  10).  Moreover,  Morkovin  (Ref.  10)  states  that 
the  three  modes  appear  to  be  uncorrelated  in  the  free  stream  of  super¬ 
sonic  tunnels;  that  is, 

lira  21  Hatf  ~  H;rr  -  0 
The  hot-wire  equation  then  reduces  to 


(e)2  . 

Aer“(f)"  .  t  \pn~{Tf)~ 

or 

Y2  - 

(*)-’  .  (7)-\-  *  (W^K,  +  n,,\)2 

where 

\  - 

^  and  V  -  T 
\cf,  \ea 

and 

Ac*77 

K 1  1  g2X 

with 

-(y-lHU2-  j) 

K 1 

1  4  (y-l)M2 

y(  1  4  n  jj\1)-(M a-l)(nt  M) 

#2 

i  4  (  0M  * 

Hot-wire  measurements  provide  the  information  needed  to  compute 
values  for  X  and  Y  in  plotting  "direct”  mode  diagrams.  If  a  value  for 
nx  is  assumed,  the  equation  can  be  graphically  solved  for  values  of 
T>  °>  and  **  using  methods  given  in  detail  by  Morkovin  (Ref.  10).  The 
reliability  of  these  procedures  is  improved  by  increasing  the  number  of 
measurements,  that  is,  the  number  of  distinct  wire  heating  currents 
for  which  data  are  obtained,  to  the  extent  that  the  additional  points  im¬ 
prove  upon  the  accuracy  of  the  fairing  curves. 


3.2  RESULTS 

Measurements  with  a  constant-current  hot-wire  anemometer  were 
made  in  the  free  stream  of  the  Tunnel  D  test  section  at  Mach  number  4 
for  three  unit  Reynolds  numbers:  0.  05,  0.  12,  and  0.  24  million  per  inch 
The  measurements  were  made  with  the  wire  positioned  on  the  tunnel 
centerline  at  the  centerline  of  the  test  section  windows,  55.  1  in.  from 
the  nozzle  throat,  and  parallel  to  the  flexible  plates  of  the  nozzle.  In¬ 
formation  regarding  the  nature  of  the  Tunnel  D  wall  boundary  layers  has 
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been  presented  by  Bell  in  Ref.  18.  From  Bell's  measurements  it  was 
determined  that  the  wall  boundary  layers  at  the  window  centerline  were 
transitional  for  Re/in.  =  0.  05  million  and  were  fully  turbulent  for 
Re /in.  =  0.  12  and  0.  24  million. 

The  hot-wire  measurements  from  the  present  investigation  have 
been  presented  in  various  forms  in  Figs.  6  through  10.  The  data  were 
obtained  using  the  input  transformer  set  on  the  higher  frequency  range 

in  40  11  w  ;  (flgnais  recorded  on  the  lower  frequency  range  (25  Hz 
to  40  kHz)  of  the  transformer  were  in  general  identical  in  rms  voltage 
uctuation  level  to  those  recorded  on  the  higher  range.  )  The  low- 

irRnT5nCy  ^'0ff  Was  set  at  1  Hz  and  the  high-frequency  cut-off  at 
lbO  kHz  The  time  constant  of  the  wire  was  determined  by  the  square- 
wave  technique  for  each  measurement  point.  Each  measurement  was 
subsequently  corrected  for  the  electronic  noise  of  the  hot-wire  instru¬ 
mentation.  The  corrected  data  were  reduced  in  accordance  with  the 
equations  presented  in  the  preceding  section. 

The  data  are  presented  in  Fig.  6  in  the  form  of  the  Kovasznay-type 
o  mode  diagram,  that  is,  e/AeT  plotted  versus  Aem/AeT.  For  each 
unit  Reynolds  number,  a  straight  line  is  the  best  simple  fairing  for  the 
ata,  and  this  is  in  accord  with  the  findings  of  Laufer  (Kef.  4)  for  free- 
stream  flow  fluctuations  in  the  working  section  of  a  supersonic  wind 
tunnel.  Data  for  wire  overheat  values  (afr  =  <RW  -  Rr)/Rr  <  0.  56)  have 
een  included  in  these  diagrams.  The  measurements  at  higher  overheat 
were  not  used  because  the  temperature  sensitivity  AeT  was  a  small  posi¬ 
tive  number  or  a  negative  number,  and,  consequently,  the  coordinates 
were  quite  large  (with  an  amplified  uncertainty)  or  negative.  This  diffi¬ 
culty  is  mathematical  in  nature  and  pertains  to  the  Kovasznay-type  of 

thlXfr18^"1,  MDat,a  at  high  overheat  values  may  be  handled  by  using 
the  direct  or  Morkovm-type  of  mode  diagram.  8 

By  following  the  arguments  of  Laufer  (Ref.  4).  the  only  simple 
fluctuating  field  which  is  consistent  with  the  straight-line  (positive 
slope)  relation  is  a  pure  sound  field  in  which  the  isentropic  relation 
between  the  fluctuating  quantities  holds.  Therefore, 


Ap  =  A _y_  AT 


and 


Au 


y~  i  t 
]  Ap 

y  p 


=  aM(y-I)(A£  +  M2M 

T0  VP  “/ 
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The  basic  equation  can  now  be  written 

37“  -  f  \r4^  -  ^  Ap  ,  M 

T  M  V1’  "  /  a-t  \y p  u  ) 


The  mean-square  voltage  is  now 

(Ap)2  Q  ’^''m  2  Au2  J  ,  Ap  ~|  r  Ap_1  ~ 

*  p "  ^ f  pH) "  *#-  ■ 

For  the  straight  line  relation  in  the  flow  diagram. 


(ApT 
A  p  'p 


aM<y_l)  "  777 


H.,p  = 


Au  Ap 


-I 


up 

so  that  rms  voltage  can  be  expressed  by 

^  'FN  1 •  (h  -  |;f; 

slope  o  the  data  fairing.  The  fluctuations  of  density  ™nd  velocity  have 

turn  ti  n  °mA  the  intercePt  and  slope  information,  and  these,  in 
turn  have  been  used  to  determine  the  magnitude  of  the  other  fluctuating 
quantities  all  of  which  are  shown  in  Fig.  7\  remembering  That 

TUP  "  T  *‘u  Seen  that  a11  of  these  fluctuation  measurements  in- 

2rs^tSK.tMM"1  wal1  boundary  layer  thickness  (d— ^ 

1  he  rms  pressure  fluctuations,  normalized  by  the  free-stream 

Rey nTlds Tifm b e T  ^  T*  ^  Fig'  »2  ^  versus  unT 
the  18  by  pTT  P°lntS  obtained  Laufer  (Ref  4)  from 

M  4  are  also  ;hUPerST1C  T™*1  °f  the  Jet  ProPulsio"  Laboratory  at 
“  are  also  shown.  For  both  sets  of  data,  the  measured  mean- 

square  values  have  been  divided  by  four,  so  that  the  data  shown  in 

»a  wa.rT,uf°XT1n-rfdiH'‘°na  **•  layer  of  one 

.  •  aufer  (Ref-  5)  Lsts  three  assumptions  which  are  required 

to  justify  this  step:  one,  that  the  four  tunnel  wall  boundary  layers 

In  Pigs.  7,  8,  and  9,  an  indication  has  been  given  of  the  possible 
effects  on  the  fluctuation  parameters  of  a  change  in  the  fairing^f  the 
mode  diagrams  A  possible  data  spread  is  shown  whieh  was  determined 
by  considering  the  maximum  and  minimum  possible  slopes  of  a  lToT 
fairing  among  the  data  for  each  unit  Reynolds  number  (Fig  G) 
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contribute  equally  to  the  sound  measured  at  a  given  point  in  the  free 
stream;  two,  that  no  correlations  exist  among  the  four  radiated  fields, 
and  three,  that  no  reflections  occur  from  opposite  walls.  It  is  noted  in 
Fig.  8  that  the  rms  pressure  fluctuations  (p/q^)  measured  for  the 
smaller  ( 12-  by  12-in. )  VKF  Tunnel  D  are  lower  than  those  for  the  JPL 
{18-  by  20-in.)  tunnel.  In  either  case,  the  fluctuation  (sound)  level  in¬ 
creased  with  increased  turbulent  tunnel  wall  boundary  layer  thickness. 

Laufer  (Ref.  5)  found  that  the  Reynolds  number  effect  on  pressure 
fluctuations  could  be  removed  by  normalizing  by  the  average  wall  shear¬ 
ing  stress.  The  measurements  from  the  present  investigation  have 
been  handled  in  this  manner  using  values  of  shearing  stress  obtained 
from  theoretical  skin  friction  values  determined  by  Van  Driest 's  second 
method  (Ref,  19)  to  normalize  the  pressure  fluctuations,  and,  as  can  be 
seen  in  Fig.  9,  the  Reynolds  number  effect  was  almost  entirely  re¬ 
moved.  Laufer's  Mach  number  4  results  (from  Fig.  12  of  Ref.  4)  have 
been  treated  in  the  same  manner,  and  these  results  are  also  shown  in 
Fig.  9. 

By  assuming  that  the  fluctuation  measurements  have  been  made  in 
the  far-field  of  the  sound  source  (that  is,  "many"  wavelengths  of  the 
field  away  from  the  source),  the  sound  waves  may  oe  considered  to  be 
plane  (Ref.  4).  The  particle  velocity  fluctuation  normal  to  the  wave 
front  (|Aun|)  is  related  to  the  pressure  variations  by  the  equation 

lAunl  Ap 

TT  M  yp 

Inasmuch  as  a  hot  wire  is  sensitive  to  velocity  variations  only  in  the 
direction  of  the  mean  flow,  it  follows  that 

Au  -  |Aun|  •  n x 

and 

Au  1  Ap 

li  M  yp  n* 

where  nx  is  the  direction  cosine  of  the  normal  to  the  wave  front  rela¬ 
tive  to  the  direction  of  the  mean  flow.  The  rms  fluctuations  can  be 
related,  with  RUp  =  -1,  by  the  equation, 

TT  1  'p 

—  =  -  —  nx 
u  M  yp 

For  the  case  of  a  stationary  Mach  wave,  nx  -  -1/M;  then  the  rms  pres¬ 
sure  fluctuation  associated  with  a  fluctuating  Mach  wave  is  determined 
from  the  equation  above  to  be 
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Ihc  equation  for  rms  voltage  with  this  substitution  becomes 


which  is  a  straight  line  passing  through  the  origin.  Inasmuch  as  the 
straight  lines  of  the  mode  diagrams  of  the  present  data  do  not  pass 
through  the  origin  (see  Fig.  13  of  Ref,  1),  it  is  concluded  that  the 
related  sound  is  not  (primarily)  radiated  by  unsteady  Mach  waves  which 
originate  from  imperfections  of  the  tunnel  nozzle  contour,  joints,  and 
similar  generators,  the  so-called  "shivering  Mach  waves"  of  Morkovin 
(Refs.  2  and  20/.  The  possibility  cannot  be  overlooked,  however,  that 
shivering  Mach  waves  may  co-exist  with  other  sources  of  sound  (see 
Section  4.  1). 


With  the  idea  of  a  fixed  disturbance  as  the  primary  source  of  the 
sound  fluctuations  ruled  out,  the  concept  of  sound  radiated  by  unsteady 
Mach  waves,  which  originate  from  a  moving  source  (Refs.  4  and  21), 
has  been  adopted, 


With  regard  to  Mach  waves  moving  relative  to  the  tunnel  walls, 
if  it  is  assumed  that  all  sound  fluctuations  are  generated  from  the  same 
source,  a  preferred  orientation  can  be  derived.  For  this  purpose, 

Laufer  (Ref.  21)  has  adapted  the  theory  of  Phillips  (Ref.  22)  for  the 
radiation  of  sound  by  a  free  shear  layer  to  the  case  of  boundary  layer 
sound  radiation.  In  essence,  a  frozen  eddy  pattern,  like  a  wavy  wall, 
moving  downstream  within  the  boundary  layer,  with  a  velocity  relative  to 
which  the  free-stream  velocity  is  supersonic,  radiates  energy  in  the  form 
of  Mach  waves.  If  the  mean  propagation  velocity  (streamwise)  of  the 
sound  source  is  designated  by  us,  the  mean  velocity  of  the  free  stream 
relative  to  the  virtual  wavy  wall  is  Gr  =  u  -  us,  and  the  corresponding 
Mach  number  is  Mr  =  Gr/a  =  (u  -  u8)/a  =  (-l/nx)  >  1.  Then,  with 
*up  = 


The  propagation  velocities  calculated  in  this  manner  for  the  pres¬ 
ent  hot-wire  measurements  are  shown  in  Fig.  10,  with  two  published 
values  determined  at  Mm  =  4  by  Laufer  (Ref.  4)  for  the  18-  by  20-in. 
tunnel  of  the  Jet  Propulsion  Laboratory.  In  comparison  with  the  data 
from  the  JPL  tunnel,  the  data  from  the  VKF  12-  by  12-in.  tunnel  exhib¬ 
ited  a  similar  variation  of  sound-source  propagation  velocity  with  chang¬ 
ing  unit  Reynolds  number  at  =  4,  but  the  source  velocity  ratio  was 
higher  for  the  smaller  tunnel  for  the  range  of  unit  Reynolds  number  in¬ 
vestigated.  The  corresponding  angle  between  the  normal  to  the  wave 
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front  and  the  direction  of  the  flow  was  between  122  and  128  deg  for  the 
three  values  of  unit  Reynolds  number  of  the  present  data.  If  this 
orientation  of  the  sound  source  is  accepted  for  the  present  measure¬ 
ments,  the  hot  wire  on  the  tunnel  centerline  would  not  sense  fluctuations 
radiated  from  the  tunnel  wall  boundary  layers  which  originated  along  the 
wall  at  less  than  7  to  7.  5  in.  upstream  of  the  window  centerline. 

It  is  seen  from  this  discussion  that  the  measurements  of  the  pres¬ 
ent  study  are  consistent  with  the  assumption  of  a  sound  source  moving 
streamwise  within  the  tunnel  wall  boundary  layer  and  radiating  energy  in 
the  form  of  Mach  waves. 

With  reference  to  the  reasoning  of  Morkovin  (Ref.  2)  that  vorticity 
and  entropy  fluctuations  in  the  test  section  have  their  origin  upstream  of 
the  nozzle  (see  Section  I),  the  following  additional  information  is  supplied 
relative  to  Tunnel  D:  (1)  An  estimate  was  made  for  Tunnel  D  of  the 
velocity  fluctuations  convected  from  the  settling  chamber  into  the  test 
section  at  Mach  number  4,  based  on  the  work  of  Ribner  and  Tucker 
(Ref.  23).  The  results  indicate  that  rms  velocity  fluctuations  are  re¬ 
duced  to  less  than  0.  004  percent  of  their  magnitude  in  the  settling 
chamber.  (2)  The  temperature  of  the  air  supplied  to  Tunnel  D  is  not 
directly  controlled  but  is  dependent  on  the  temperature  to  which  the  air 
is  heated  in  the  high  pressure  air  storage  reservoir  and  on  the  pressure 
drop  across  the  control  valves  which  regulate  stagnation  pressure.  Per¬ 
forated  baffle  plates  in  the  upstream  diffuser  between  the  valves  and  the 
settling  chamber  (Fig.  lc)  cause  the  flow  to  mix,  and  the  screens  in  the 
settling  chamber  further  promote  the  decay  of  any  temperature  spots 
(entropy  fluctuations)  which  are  convected  into  the  supply  section.  It  is, 
therefore,  concluded  that  vorticity  and  entropy  fluctuations  in  the  test 
section  of  Tunnel  D  are  negligibly  small  and  that  the  assumption  of  a  pure 
sound  field  is  well  founded. 


SECTION  IV 

DISCUSSION  OF  HOT-WIRE  POWER  SPECTRA 


4.1  SPECTRAL  ANALYSIS  OF  HOT  WIRE  DATA 

Further  sturcture  information  that  characterizes  the  fluctuations 
in  the  flow  field  is  obtained  through  the  energy  spectrum  density  of  the 
hot-wire  output.  In  this  type  of  data  analysis,  the  distribution  of  the 
kinetic  energy  of  the  fluctuations  with  frequency  or  wave  number  is  given. 
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ceding  section  fw  b  dd?d  h®  separating  techniques  of  the  pre- 

-  ttZSritfz a—  - 
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theses  denote  functional  relationships  X As  hrfbSCr*^'  and  the  paren_ 
based  on  an  assumption  of  a  sin./i*  h  •  bcfore,  this  equation  is 
ence  of  all  three  modes  and  th  "n  d°Tf?lt  souad  source.  In  the  pres- 

require  at  least  six  different  overheat's aTaTman^ fCallbrati0n  would 
possible,  For  the  case  nf  th*  ™  *  t  S  m  ny  frequencies  as 

the  dominance  of  the  sound tS ^sent  "^surements,  the  data  indicate 
and  the  equation  reduces  to  ^  “  dlSCUSS6d  in  the  prcced-S  sections, 

K"(f-A»)  -  Ar2(A;,0K^0 
for  a  single  sound  source. 

Ideally,  the  energy  spectral  density  is  obtained 


as 
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assumption  that  the  fluctuation  «=•  i  .mitely  lonS  time  (t).  Under  the 
tom  of' „„ro“  2»Z7JlTU°",«y  -  “««“'■  •  W- 

(Af)  can  be  used  to  obtain  the  so  /  V®*  (t)  and  a  finitc  bandwidth 

car  error  „  ■«**- 


errur 


K(f,  \‘)  ’ 
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where  s  is  the  standard  deviation  and  E(f  A'  )  is  th„ 

termined  spectrum.  For  the  case  JhLV  experimentally  de- 

Af  =  50  Hz  and  t  =  4  sec  yielding  a  stnt  “Cafuroments  reported  herein, 
cent.  '  y  cldlnS  a  statistical  uncertainty  of  ~±7.  0  per- 
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Further  complicating  the  experimental  determination  of  the  energy 
density  function  is  the  matching  of  the  anemometer  compensation  ampli¬ 
fier  to  the  wire  frequency  roll-off  so  that  the  wire  has  an  effective  flat 
frequency  response.  This  is  accomplished  through  the  time  constant 
setting  in  the  amplifier.  If  the  time  constant  is  properly  set,  the  fre¬ 
quency  response  is  as  indicated  in  Fig.  11a.  However,  as  was 
mentioned  in  the  preceding  section,  because  of  experimental  difficulties, 
a  degree  of  uncertainty  in  the  selection  of  the  time  constant  was  inherent 
in  the  experiment.  The  effect  of  improper  compensation  on  the  wire  re¬ 
sponse  is  shown  in  Figs,  lib  and  c.  As  a  result  of  this  effect  some  dis¬ 
tortion  of  the  measured  spectra  can  be  anticipated. 

In  interpreting  the  measured  spectra  when  the  fluctuations  are 
considered  to  be  sound -induced,  the  relative  propagation  of  the  sound 
mode  introduces  a  frequency  distortion  when  sensed  by  a  stationary  hot 
wire.  This  is  illustrated  by  considering  a  plane  wave  with  the  wave 
vector  inclined  to  the  mean  flow  direction  by  an  angle  9.  The  frequency 
a  +  u  cos  Q  a 

is  given  by  f  = - - -  =  (1  +  M  cos  6),  and  when  cos  0  8  -1/M,  a 

Mach  wave,  the  frequency  indicated  by  the  wire  is  zero.  Consequently, 
it  can  be  inferred  that  frequencies  near  zero  represent  Morkovin's 
"shivering  Mach  waves"  (Refs.  2  and  20).  For  this  special  case,  the 
corresponding  mode  diagrams  would  pass  through  the  origin,  as  was 
pointed  out  in  the  previous  section,  and  the  equation  for  the  sensitivity 
coefficient  Ae^  (see  Section  3.  1)  would  use  the  substitution  nx  -  -  1/M. 
However,  a  non-zero  value  of  the  intercept  indicates  a  traveling  sound 
source  with  a  dominant  orientation  other  than  that  of  a  Mach  wave.  In 
this  case,  the  applicable  direction  cosine  is  that  calculated  from  the 
average  source  velocity.  For  computational  purposes,  the  model 
allows  only  one  sound  source  (otherwise  a  second  sound  sensitivity 
coefficient  must  be  included).  It  must  be  pointed  out  again  that  this 
does  not  rule  out  the  possibility  that  the  real  physical  situation  may 
include  both  "shivering  Mach  waves"  and  a  traveling  sound  source. 

A  further  complicating  factor  in  interpreting  the  spectra  occurs 
because  of  the  directional  sensitivity  of  the  wire  through  the  overheating 
parameter  (A^,).  Inspection  of  the  sound  sensitivity  coefficient  Ae^ 
reveals  (Ref.  11)  that  a  value  of  A'w  can  be  found  for  which  Aen 
(Table  I)  is  independent  of  sound  direction.  In  reporting  the  results 
herein,  the  data  have  been  grouped  according  to  the  alternate  overheat 
(a^,)  and  reported  for  the  various  unit  Reynolds  numbers. 

A  more  nearly  precise  comparison  of  spectra  could  have  been 
made  with  a  constant  sensitivity  coefficient  ratio  (Aem/Aex).  With  this 
coefficient  ratio  maintained  constant,  the  sensitivity  of  Ae^  to  the  sound 
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wave  direction  would  be  constant.  This  technique  was  not  utilized  in 
and  the?  Stfdy  bt’CaUSC  thc  rcnuired  calculations  are  quitc  lengthy 
“,uW,iT„^d«aPr06r:m  *“  bCine  dCVolo|*■,  with' the 

w  „  A"  exemplary  variation  of  Aem/AcT  with  certain  nominal  values 
of  overheat  and  with  unit  Reynolds  number  is  shown  below: 


Nominal  Overheat,  aw  = 

=  0.  05 

0.4 

0.  5 

Re/in.  =  0.  05  x  10G 

0.  134 

1.  57 

2.47 

=  0.  12  x  10G 

0.  120 

1.37 

2.  71 

=  0. 24  x  10G 

0.  104 

1.  50 

2.  70 

in  thnIh„C  f^c.e'strcam  measurements  indicate  a  useful  energy  content 
m  the  signal  to  a  frequency  of  about  80  kHz.  At  this  point,  the  sicnal 
approached  ,„c  noKo  level  of  the  h«,-»lrc  a.  shTnt 

tern  is  h  !Wo^eovcJ’*  *f  thc  uPPer  frequency  limit  of  the  hot-wire  sys- 
0  707  (3  dW  f  i°qTnCy  at  which  tho  si6nal  ^  attentuated  to 

“at  ”,  (= 5oS  nf  fh  amP  l  f  at  ZCr°  frcciucnc>‘  »ho  "ceiling -to -floor" 
unnnr  -  5°°  f  *5?  ampllflcr  imposes  a  lower  bound  for  the  various 
pper  ii  equency  limits  of  this  study  (wherein  the  time  constants  had  an 

ZTZlT!  m"C’  *  ~80  kHz-  ™*  ” m“t  £  n 

iTpie  12  'ft!8  of  ‘he  roll-off  of  the  hot-wire  signal 

th  k6;  12 '  thC  systom  limitation  is  based  on  the  3-db  point,  whereas 

f  =C  lOWfe”  S1Snal  ^  80  kHZ  WaS  dOW"  by  24  db  from  its  lcvcl  at 

RianaiInSP!Ctl*°n  °i  thC  anaiysis  &ot8  indicated  that  the  most  significant 
gnal  fluctuation  frequencies  were  below  20  kHz,  while  above  this  fre¬ 
quency  all  spectra  appear  to  decay  in  an  orderly  manner  for  this  mason 

Xs'uEr  °f  ”*ar'  WOW  M  kHZ  *S  <"  'he  shooed’. 

that  ThC  Te!hu  d  °f  spectral  data  Presentation  selected  for  Fig  13  is 
that  suggested  by  Liepmann,  Laufer,  and  Liepmann  (Ref.  25)  in  which 

[*<o|  df  w  [Tr  J'W Odf 


where 


r~i2 

is  thr  sound  cnerg\  density  at  f 
[ff]"  is  Uie  lotul  sound  energ) 
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^»(f)  in  the  normalized  sound  energy  density  distribution  function 

satisfying  the  condition  /°°Uf)df  *  1 
o 

The  particular  advantage  of  this  form  lies  ir,  'he  elimination  of  the  cali¬ 
bration  coefficients,  since 


Mf) 


-  0(f) 


The  normalized  sound  energy  density  is  obtained  from  the  energy 
density  distribution  by 


-  M2G(f) 


where  ?  is  obtained  as  indicated  in  Section  3.1.  This  particular  data 
presentation  technique  is  less  clear  whenever  there  is  an  uncertainty  of 
the  direction  cosine  nx  of  the  sound  wave. 


As  a  final  note  of  caution  in  interpreting  the  data,  the  above  tech¬ 
niques  have  specifically  ruled  out  the  existence  of  other  modes;  that  is, 
entropy,  vorticity,  or  another  sound  source.  The  conclusion  that  they 
contribute  only  negligibly  to  the  turbulence  energy  is  reasonably  justi¬ 
fied,  but  on  the  other  hand,  they  may  contribute  significantly  as  an  effect 
on  other  phenomena. 


4.2  RESULTS 

Power  spectral  density  analyses  of  the  hot-wire  signal  were 
obtained  for  three  nominal  wire  overheats,  a^,  =  0.  05,  0.  4,  and  0.  5, 
at  each  of  the  three  unit  Reynolds  numbers  of  the  present  study,  and  the 
results  are  presented  in  various  forms  in  Figs.  13  through  15.  An  un¬ 
certainty  of  ±1  db  is  estimated  for  the  power  spectral  density  traces 
(which  are  not  shown),  and  the  corresponding  uncertainty  interval  is 
indicated  in  each  figure  for  several  signal  frequencies  for  the  data  ob¬ 
tained  at  Re/in.  =  0.  12  x  106;  these  intervals  are  typical  of  all  three 
values  of  Re/ in. 

The  normalized  energy  density  distribution  for  the  nominal  over¬ 
heats  is  presented  in  Fig.  13  for  frequencies  f  <  20  kHz.  The  lack  of 
similarity  of  the  spectra  for  a  given  unit  Reynolds  number  and  different 
overheats  is  a  result  of  the  dependence  of  the  sound  sensitivity  on  over¬ 
heat  (Ref.  16).  The  differences  in  shape  at  each  overheat  warrant  dis¬ 
cussion  of  each  set  separately.  The  low  overheat  <aw  =  0.  05)  data  plot 
(Fig.  13a)  exhibits  sharp  peaks  at  frequencies  less  than  1  kHz.  The 
maximum  peak  is  at  the  lowest  unit  Reynolds  number,  the  lowest  peak 
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is  at  the  intermediate  unit  Reynolds  number,  and  the  highest  unit 
Reynolds  number  has  the  intermediate  peak,  thus  indicating  the  follow¬ 
ing  trend:  as  the  unit  Reynolds  number  increased,  the  amplitude  of  the 
low  frequency  portion  of  the  spectrum  (f  <  1  kHz)  at  first  decreased  and 
then  increased.  A  similar  trend  is  exhibited  at  the  higher  overheat, 
aw  =  0.  4,  where  the  signal  amplitude  is  much  larger.  Inspection  of  the 
plots  for  all  three  overheats  in  Fig.  13  reveals  that  the  amplitude  at 
any  given  frequency  f  >  1  kHz  was  essentially  the  same,  within  experi¬ 
mental  error,  for  both  the  intermediate  and  highest  unit  Reynolds  num¬ 
bers.  When  the  normalized  sound  energy  density3  is  plotted,  as  in 

^-4,  it  is  seen  that  in  general  the  normalized  soind  energy  density 
decreased  with  increasing  unit  Reynolds  number  in  the  low  fre¬ 
quency  range.  It  should  be  mentioned  that  this  low  frequency  effect  is 
much  more  pronounced  than  is  indicated  by  the  mean  square  fluctuation 
quantities.  If  the  data  for  Re/in.  =  0.  05  x  10^  are  taken  as  a  reference, 
the  percentage  change  of  the  peak  values  for  the  remaining  unit  Reynolds 
numbers  varies  from  60  to  89  percent  for  the  three  overheats.  This  is 
in  contrast  to  a  maximum  change  of  50  percent  in  the  mean  square 
fluctuation,  (?)3. 

Laufer  (Ref.  5)  and  Laufer,  Ffowcs  Williams,  and  Childress 
(Ref.  26)  have  found  that  the  use  of  an  integral  scale  to  normalize 
’’far  field"  energy  density  spectra  at  Mach  numbers  2.  0  and  4.  5  pro¬ 
duced  very  similar  distributions  for  the  two  Mach  numbers.  It  can  be 
concluded  that  in  addition  to  the  mean-square  fluctuation  which  char¬ 
acterizes  the  flow,  an  additional  parameter,  such  as  a  length  scale,  is 
needed  to  characterize  the  spectral  shape.  In  lieu  of  the  integral  scale 
as  developed  in  Ref.  5,  the  thickness  (6)  of  the  tunnel  wall  boundary 
layer  was  chosen  for  normalizing  the  present  data  for  a  similar  presen¬ 
tation  (Fig.  15),  The  values  used  for  <5  were  taken  from  the  measure¬ 
ments  made  by  Bell  (Ref.  18)  at  the  window  centerline  station  of 
Tunnel  D. 

The  normalized  energy  spectra  for  the  three  Reynolds  numbers  of 
the  present  study  at  an  overheat  of  a^  =0.4  have  been  presented  in 
Fig.  15  in  the  form  u^Gff)/ 6  plotted  against  f  6/ uxJ  for  frequencies 
250  Hz  <  f  <  80  kHz.  Data  for  the  two  higher  values  of  unit  Reynolds 
number  are  in  excellent  agreement  for  f  6/u.  >  0.  035  (that  is  for 
f  >  700  Hz).  Moreover,  at  lower  frequencies  where  the  time  constant 
of  the  analysis  equipment  is  relatively  long,  the  uncertainty  of  the  re- 
sults  is  probably  greater  than  indicated;  in  which  event,  Mie  separation 

3It  has  been  assumed  that  the  spectra  arc  sound  dominaieu  end 
result  from  a  single  source  with  an  average  clirectio;  eoeine  deter¬ 
mined  from  the  mode  diagram  analysis. 
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of  the  two  curves  may  lack  significance.  With  a  decrease  in  unit 
Reynolds  number  to  0.  05  x  106  per  inch,  more  energy  becomes  concen¬ 
trated  at  the  lower  frequencies,  while  for  intermediate  frequencies 
(3  kHz  <  f  <  20  kHz)  the  differences  among  the  three  curves  was  not  sig¬ 
nificant.  For  higher  frequencies  (f  >  20  kHz),  the  level  of  the  hot-wire 
signal  for  the  lowest  unit  Reynolds  number  (0.  05  x  106  per  inch) 
approached  the  level  of  the  electronic  noise  of  the  hot-wire  equipment 
(within  3  db)  and,  therefore,  was  considered  to  be  no  longer  meaning¬ 
ful  in  this  frequency  range. 

The  explanation  offered  by  Laufer,  Ffowcs  Williams,  and 
Childress  (R.ef.  26)  for  the  behavior  of  pressure  spectra  in  the  far 
field  was  obtained  through  the  "eddy  Mach  wave"  concept.  In  this 
model,  the  eddy  scales  in  the  turbulent  boundary  layer  determine  the 
scales  (wavelengths)  of  the  radiated  sound.  For  an  increasing  unit 
Reynolds  number,  the  wall  boundary  layer  thickness  decreases  with  the 
result  that  the  energy- containing  eddies  become  smaller,  and  this 
forces  a  shift  of  the  energy  density  spectra  to  higher  frequencies. 
Further  investigation  should  be  made  of  the  power  spectra  of  the  radi¬ 
ated  sound  mode  so  that  the  disturbances  can  be  more  nearly  completely 
characterized.  J 


SECTION  V 

DISCUSSION  OF  SURFACE  PRESSURE  FLUCTUATIONS 
ON  A  FLAT  PLATE 


In  order  to  relate  the  tunnel  wall  radiated  sound  measured  in  the 
free  stream  by  the  hot-wire  technique  to  pressure  fluctuations  mea¬ 
sured  on  a  model  surface,  a  flat  plate  instrumented  with  a  flush-mounted 
microphone  was  installed  in  the  test  section  of  Tunnel  D  at  the  conclusion 
of  the  hot-wire  measurements.  The  8 -in. -long,  5-in. -wide  flat  plate 
(Fig.  3)  used  by  Pate  and  Schueler  (Ref.  6)  was  mounted  with  the  flat 
surface  located  1.2  in.  above  the  tunnel  centerline  (for  <j>  =  0)  and  with 
the  microphone  positioned  axially  at  the  window  centerline. 

Data  were  obtained  at  the  following  nominal  tunnel  conditions: 


psia  Re/in.  x  10"c 


3 

0.  025 

6 

0.  050 

15 

0.  12 

30 

0.  24 

60 

0.  48 
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The  results  normalized  by  the  free -stream  dynamic  pressure  are 
presented  in  Fig.  1G.  Similar  measurements  obtained  with  the  flat  plate 
rotated  counterclockwise  about  the  tunnel  centerline  (<J  =  -90  deg)  have 
been  included  also.  The  hot-wire  anemometry  results  from  the  present 
study,  apportioned  for  the  radiation  of  a  single  tunnel  wall  (Fig.  8),  are 
shown  again  in  Fig.  1G,  for  comparison  purposes.  Regarding  the  micro¬ 
phone  data,  it  is  noted  that  the  measurements  at  Rc/in.  =  0.  12  x  10^, 
0.24x10^,  and  0.48  x  10  ®  are  effectively  equal.  For  Rc/in.  =  0.05  x  10  ,  the 
pressure  fluctuations  were  generally  lower;  whereas  at  Re/in.  =  0.025  x  10G, 
the  measurements  were  an  order  of  magnitude  lower.  The  data  obtained 
with  the  microphone  oriented  toward  the  tunnel  sidewall  (<t>  =  -90  deg) 
were  in  general  somewhat  lower  than  those  recorded  with  the  micro¬ 
phone  oriented  toward  the  flexible  plate  (<£  =0),  but,  by  noting  the  spread 
in  the  measurements  for  =  0,  these  differences  may  lack  significance. 

In  the  calculations,  no  attempt  was  made  to  correct  the  data  for  the 
attenuation  of  microphone  response  which  would  be  expected  because  of 
the  finite  size  of  the  microphone  diaphragm  (0.  25-in. -diam)  compared 
with  the  short  wave  lengths  contained  among  the  fluctuations  (Refs.  27 
and  28) . 

The  magnitude  of  free -stream  pressure  fluctuations  at 
Re/in.  =  0.  025  x  10G  would  be  of  interest  for  a  comparison  with  the 
microphone  data;  unfortunately,  the  hot-wire  measurements  were  not 
extended  below  a  value  of  Re/in.  of  0.  05  million.  However,  Vrebalovich 
has  reported  (Ref.  29)  that  free-stream  fluctuation  levels  were  found  to 
be  smaller  when  the  tunnel  wall  boundary  layer  was  laminar  than  when 
it  was  turbulent.  Therefore,  with  decreasing  unit  Reynolds  number,  a 
fall-off  in  free-stream  pressure  fluctuation  level  somewhere  below  the 
value  for  Re/in.  of  0.  05  million  would  be  expected  from  hot-wire  mea¬ 
surements  in  Tunnel  D. 

The  Tunnel  D  wall  boundary  layer  thickness  measurements  by 
Bell  (Ref.  18)  indicate  that  the  flexible  plate  boundary  layer  was  tran¬ 
sitional  at  the  window  centerline  station  for  values  of  unit  Reynolds  num¬ 
ber  between  0.  04  and  0.  08  million  per  in.  (By  using  the  orientation  of 
the  sound  source  derived  in  Section  3.  2  and  noting  again  that  the  flat- 
plate  surface  was  located  1.2  in,  off  tunnel  centerline,  it  follows  that 
the  microphone  would  not  receive  sound  radiation  which  originated 
along  the  tunnel  wall  at  less  than  5.  5  to  G  in,  upstream  of  the  window 
centerline.  )  It  is,  therefore,  concluded  that  the  present  microphone 
data  for  Re/in.  =  0.025  and  0.05  million  were  associated  with  sound 
radiated  by  laminar  and  transitional  tunnel  wall  boundary  layers,  re¬ 
spectively. 
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Moreover,  if  one  considers  the  findings  of  Pate  and  Brown 
(Ref.  31)  and  Pate  and  Schueler  (Figs.  5  and  6  of  Ref.  6),  it  could  be 
inferred  that  the  maximum  sound  radiation  would  be  associated  with  the 
completion  of  transition;  in  which  event,  the  free-stream  fluctuations 
probably  attained  a  higher  level  than  any  actually  measured  in  the  pres¬ 
ent  investigation,  that  is,  between  present  values  of  unit  Reynolds  num¬ 
ber.  This  reasoning  led  to  the  fairing  included  in  Fig.  16,  although 
the  level  and  precise  unit  Reynolds  number  of  the  maximum  are  neces- 
sarily  arbitrary,  here. 

A  comparison  of  the  two  sets  of  measurements  shown  in  Fig.  16 
reveals  that  the  pressure  fluctuations  determined  by  the  microphone 
under  the  boundary  layer  of  the  flat  plate  were  20  to  45  times  larger 
than  those  determined  in  the  free  stream  by  the  hot  wire.  This  magni¬ 
tude  of  amplification  of  pressure  fluctuations  by  the  flat  plate  boundary 
layer  is  consistent  with  the  recent  findings  of  Kendall  (Ref.  31)  at  a 
Mach  number  of  4.  5  in  the  JPL  18-  by  20-in.  tunnel.  The  implication 
is  that  the  laminar  boundary  layer  over  its  3.  5 -in.  length  is  responsible 
for  at  least  a  20-fold  amplification  or  accumulation  of  free-stream  pres¬ 
sure  fluctuations.  It  should  be  noted  again  that  the  microphone  measure¬ 
ments  have  been  presented  without  any  corrections  for  attenuation  of 
response  resulting  from  the  finite  size  of  the  sensor  and  are,  therefore, 
conservative. 

A  power  spectral  density  analysis  of  the  microphone  output  re¬ 
corded  at  Re/in.  =  0.  24  x  106  is  presented  in  Fig.  17  which  is  representa¬ 
tive,  except  for  signal  level,  of  all  analyses  of  microphone  signals  re¬ 
corded  in  the  present  investigation.  The  output  signal  was  analyzed 
using  a  bandwidth  of  50  Hz  for  the  frequency  range  shown  in  the  figure. 

The  trace  is  characterized  by  several  outstanding  peaks  which,  in 
general,  were  found  in  all  of  the  analyses.  The  most  prominent  of  the 
peaks,  at  f  =  45  kHz,  is  the  result  of  the  microphone  resonance  pre¬ 
viously  discussed.  The  sources  of  the  remaining  peaks  are  unknown, 
but  signals  from  the  hot  wire  in  the  free  stream  for  the  same  tunnel  con¬ 
dition  do  not  exhibit  these  spikes.  Whether  they  are  related  to  the  micro¬ 
phone,  to  the  amplification  by  the  flat  plate  boundary  layer,  or  to  some 
other  cause  should  be  investigated. 

The  power  spectra  for  a  set  of  microphone  data  points  which 
spanned  the  range  of  unit  Reynolds  numbers  from  Re/in.  =  0.  025  x  106 
to  0.48  x  10  are  compared  in  Fig.  18.  The  similarity  of  the  shapes  of 
the  microphone  energy  density  spectra,  especially  for  Re/in.  >  0.  05  x  10^, 
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was  not  expected.  This  peculiarity  also  should  be  investigated  to  deter¬ 
mine  if  it  is  an  actual  physical  phenomenon  or  was  a  measurement  idio¬ 
syncrasy.  This  presentation  shows  that,  with  increasing  unit  Reynolds 
number,  the  level  of  the  signal  at  each  frequency  rose  to  an  essentially 
constant  value  associated  with  the  tunnel  conditions  (Re/ in.  =  0.  12  x  10® 
and  0.  24  x  10®)  for  which  the  tunnel  wall  boundary  layers  were  turbulent. 
The  trace  presented  for  Re/in.  =  0.  48  x  10®,  however,  reveals  a  trend 
not  obvious  among  the  rms  pressure  fluctuation  data  (Fig.  16),  namely 
another  rise  in  signal  level  with  the  tunnel  wall  boundary  layer  already 
fully  turbulent.  An  examination  of  Schlieren  photographs  made  of  the 
flow  over  the  flat  plate  as  the  microphone  data  were  being  recorded 
revealed  that  the  flat-plate  boundary  layer  was  laminar  at  the  micro¬ 
phone  station  for  Re/in.  ^  0.  24  x  10®  and  had  become  transitional  for 
Re/in.  =  0.48  x  10®.  It  is  concluded  from  this  evidence  that  the  increase 
in  the  level  of  the  power  spectra  shown  for  the  latter  unit  Reynolds  num¬ 
ber  was  a  result  of  the  change  in  the  nature  of  the  flat-plate  boundary 
layer  above  the  microphone. 

The  lack  of  an  analytic  theory  and  the  widely  different  measure¬ 
ment  characteristics  of  the  hot  wire  and  the  microphone  leave  the 
nature  of  the  interactions  of  free-stream  and  flat-plate  boundary  layer 
fluctuations  unresolved  in  this  investigation.  An  interesting  extension  of 
this  study  could  be  performed  with  a  hot  wire  immersed  in  a  flat-plate 
boundary  layer  that  has  a  flush-mounted  microphone  so  that  both  sensors 
experience  the  same  flow  field. 


SECTION  VI 

CONCLUDING  REMARKS 


Hot-wire  anemometry  techniques  were  successfully  employed  to 
determine  free-stream  flow  fluctuations  at  a  Mach  number  of  Mffl  =  4  in 
the  test  section  of  the  12-  by  12 -in.  supersonic  tunnel  (D).  The  results 
are  quite  similar  to  findings  at  =  4  in  another  tunnel  which  have  been 
published  by  Laufer.  It  was  found  that  each  of  the  various  fluctuations 
deduced  from  the  hot-wire  measurements  increased  with  the  thickness 
of  the  turbulent  tunnel  wall  boundary  layers.  The  Reynolds  number 
effect  on  the  pressure  fluctuations  was  removed  by  normalizing  by  the 
average  wall  shearing  stress.  The  free-stream  measurements  of  the 
present  study  have  been  examined  using  the  assumption  of  a  sound  source 
traveling  in  the  tunnel  wall  boundary  layer  at  a  velocity  relative  to  which 
the  free-stream  velocity  is  supersonic.  The  data  exhibited  an  increase  of 
the  propagation  velocity  of  the  sound  source  with  increasing  unit  Reynolds 
number. 
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Power  spectral  density  analyses  of  che  hot-wire  output  indicate 
that,  in  addition  to  the  mean-square  fluctuation  characterizing  the  flow, 
an  additional  parameter  such  as  a  length  scale  is  useful  to  characterize 
the  spectral  shape. 

Supplementary  measurements  in  the  test  section,  for  the  same 
testing  conditions  as  for  the  hot  wire,  using  a  microphone  mounted 
flush  with  the  surface  of  a  flat  plate,  revealed  pressure  fluctuation 
levels  larger  than  those  of  the  free  stream  by  a  factor  of  at  least  20. 
Further  inquiry  into  the  nature  of  the  interactions  between  the  free- 
stream  fluctuations  and  the  flat-plate  boundary  layer  is  needed  to  ex¬ 
plain  this  phenomenon. 
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Fig.  1  Tunnel  D,  a  12-  x  12-in.  Supersonic  Wind  Tunnel 


Preceding  page  blank 


35 


AEDC  TR-7M43 


b  *  g 
ft  g  o 
l  ci  o 


1  TO  CJ 

J3  H 

bfi 

D 

•h  a> 

k 

33  hfi  3  w 

TO 

(A  0) 

E  k 

(ft  > 

0  0 

0)  rH 

k  HJ 
Ph  CO 

££ 

a > 

■M  43 
O  +■>  /-^ 
O  >> 
CO  k  ^ 
O  3 
<D  <H  H-> 
W  CO 
k  T3 
<d  a>  «h 
>  >  c 

U10Q) 

c  g  1)] 
cq  a>  a» 


k 

0) 

a  c n 

o  x 

y 

>»  TO 

rH  1-3 


lO 
N  CO 
rH  O 
O  O 

O  O 


O  O 
CO 


fl> 

u 

3 

(0 

09 

<D  (fl 

rH 

O  TO 
■H  > 

+J 

TO  4H 

B  0) 
0  -H 

«H  rH 

3  0> 
<  Pi 


m 

c 

CD 

a> 

k 

43  y 

V  CO 
•H 

»  &C 

c 

k  -H 
QJ  Q« 

43  E 
E  TO 
TO  Q 
43 

u  a> 

y 

tJD  C 

c  <d 

*H  rH 
rH  3 
rH  43 
•H  k 


6 


3G 


Perspective  Drawing 
Fig.  1  Concluded 
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a.  Side  View 


b.  Top  View 
Fig.  2  Hot-Wire  Probe 
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Soft-Solder 
Connection 
(Typical) 


Wire  of  Platinum 
Rhodium  Alloy 
(90/10  Percent) 


Wire  Diameter 


0.0001  in 


c.  Shadowgram  of  a  Typical  Wire  Attachment 
Fig.  2  Concluded 
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Fig.  3  Flat  Plate  Microphone  Model 
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Recovery  Factor,  t}  -  Tr/T0  Nusselt  Number,  Nu^  =  hd/k 
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a.  Nusselt  Number 


b.  Temperature  Recovery  Factor 

Fig.  5  Hot-Wire  Heat-Loss  Variation  with  Unit  Reynolds  Number,  M  =  4.0 
d  *  10  4  in.,  T0  =  530°  R 
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Bars  in  Figs.  8  and  9  indicate  data 
spread  possible  with  alternate  fairings 
of  mode  diagrams  for  each  unit  Reynolds 
number  (Fig.  6).  y 


Fig.  8  Variation  of  RMS  Pressure  Fluctuations  (Normalized  by 
Dynamic  Pressure)  with  Unit  Reynolds  Number 


Fig.  9  Variation  of  RMS  Pressure  Fluctuations  (Normalized  by 
Wall  Shearing  Strati)  with  Unit  Reynolds  Number 
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Unit  Reynolds  Number,  Re/in. 

Fig.  10  Variation  of  Sound-Source  Velocity  with  Unit  Reynolds  Number 
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Frequency,  f.  Hz  Frequency,  f,  Hz  Frequency,  f,  Hz 

Proper  Time  Constant  Setting  b.  Time  Constant  too  High  c.  Time  Constant  too  Low 

Fig.  11  Effects  of  Time  Constant  Setting  on  Anemometer  Response 
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Frequency,  f,  Hz 

Fig.  12  Typical  Hot-Wire  Response  and  Electronic  Noise  at  High  Frequencies 
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Fig.  13  Continued 
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Fig.  14  Normalized  Sound  Energy  Density  Distribution 
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Fig.  15  Hot-Wire  Signal  Spectra  for  Wire  Overheat  of  a’w  =  0.4 
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Pressure  Fluctuation 


AEDC  TR  71  143 


Unit  Reynolds  Number,  Re/in. 

Fig.  16  Comparison  of  Flat  Plate  Surface  and  Free  Stream  Pressure 
Fluctuations 
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Microphone  Resonance 
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Fig.  17  Power  Spectral  Density  Analysis  of  Microphone  Output 
Recorded  at  Re/in.  *  0.24  x  10® 
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Fig.  18  Comparison  of  Power  Spectra  of  Microphone  Output  for 
Various  Unit  Reynolds  Numbers 
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